Design considerations for table-top, laser-based VUV and X-ray free electron lasers 
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A recent breakthrough in laser-plasma accelerators, based upon ultrashort high-intensity lasers, 
demonstrated the generation of quasi-monoenergetic GeV-electrons. With future Petawatt lasers 
ultra-high beam currents of ~ 100 kA in ~ 10 fs can be expected, allowing for drastic reduction in 
the undulator length of free-electron-lasers (FELs). We present a discussion of the key aspects of 
a table-top FEL design, including energy loss and chirps induced by space-charge and wakefields. 
These effects become important for an optimized table-top FEL operation. A first proof-of-principle 
VUV case is considered as well as a table-top X-ray-FEL which may open a brilliant light source 
also for new ways in clinical diagnostics. 

PACS numbers: 41. 60.0,52. 38.Kd 



INTRODUCTION 

The pursuit of table-top FELs combines two rapidly de- 
veloping fields: laser-plasma accelerators, where the gen- 
eration of intense quasi-monoenergetic electron bunches 
up to the GeV range has been achieved 0, [3, S 0] j an d 
on the other hand large-scale X-ray free-electron lasers 
(XFELs) that are expected to deliver photon beams with 
unprecedented peak brilliance. A prominent application 
of such FEL pulses is single molecule imaging [5| . The 
proposed laser-plasma accelerator-based FELs would not 
only allow a greater availability due to their smaller size 
and costs, but also offer new features, such as pulses syn- 
chronized to the phase-controlled few-cycle driver laser 
Q for pump-probe experiments. Moreover, the X-ray 
energy of a table-top XFEL can be as large as required 
for medical diagnostics (above 20 keV jjjj]). 

The mechanism for the generation of intense (nC 
charge) quasi-monoenergetic electron pulses by laser- 
plasma accelerators requires an ultrashort, high- intensity 
laser-pulse with a length shorter than the plasma wave- 
length (on the //m-scale corresponding to gas densities 
of 10 19 cm -3 ). Due to the ponderomotive force, plasma 
electrons are blown out transversely, leaving an electron- 
free zone - the bubble - behind the laser pulse Q. These 
electrons return to the axis after half a plasma oscillation, 
thus determining the size of the bubble in the order of the 
plasma wavelength. Typically about 10 9 ...10 10 electrons 
are captured into the bubble, as found both experimen- 
tally and from scaling laws of relevant bubble param- 
eters within the framework of similarity arguments Q. 
Due to the inertial positive ion background, these elec- 
trons experience a strong electrical field gradient of up to 
TV/m. Particle- in-cell (PIC) simulations [l(| show that 



the bubble electrons form a stem that is geometrically 
considerably smaller than the bubble as seen in Fig. 1. 
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FIG. 1: (Color online) Snapshot of electron density (in units 
of 10 20 cm -3 ) from PIC simulation. The geometrical size 
of the electron- free cavity ("bubble") behind the laser spot 
corresponds to the plasma period, which is about 8 /J,m in 
this case (gas density 1.8 • 10 19 cm -3 ). The stem of the high- 
energy electrons is much shorter than the plasma period. 

For these dimensions, beam currents of the order of 
~ 100 kA (i.e. ~ 1 nC charge within ~ 10 fs) can be 
reached without the need for any bunch compressor. 
For preparation of the required gas densities, either su- 
personic gas-jets or capillaries [ll| are used, the latter 
consisting of a ~ 300 /im-thin gas-filled channel with a 
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parabolic radial ion density profile generated by a second 
laser pulse or discharge. The advantage of the capillary is 
that the laser can be guided beyond the Rayleigh length 
allowing longer acceleration distances and hence higher 
electron energies. The very recent experiments by Lee- 
mans et al. [4|, utilizing a 40— TW laser pulse of 38 fs 
duration and a gas density of 4.3 • fO 18 cm~ 3 , clearly 
show that 1 GeV electron beams can now be produced 
with capillaries. However, the measured charge of 30 pC 
is significantly below our design goal of 1 nC. In order 
to further improve the resulting current, we propose the 
scheme described below. The number of bubble electrons 
can be increased when using plasmas with higher density 
such as in Fig. 1, because the feeding process is more ef- 
ficient if more plasma electrons are present. An increased 
gas density requires shorter laser-pulses (sub-10-fs), such 
that the entire laser pulse fits into one plasma period. 
In this case, the entire laser pulse energy can be used, 
while longer pulses lose energy during self-shortening for 
entering the bubble regime. Therefore, in contrast to 
all previous laser-acceleration experiments, we propose 
to use a different driver laser, namely a ~ 5 fs-Petawatt 
laser (like the Petawatt-Field-Synthesizer PFS at MPQ 
[II]). Own PIC simulations show that such ultrashort 
lasers can capture more than 1 nC charge inside the bub- 
ble. The smaller plasma period leads to smaller bubbles, 
hence shorter bubble stems and thus higher currents (1 
nC in 10 fs). However, the expected final energy in the 
bubble regime is, according to the scaling laws, lower 
for shorter laser pulses. To overcome this limit, we also 
propose to use a capillary and suggest the following sce- 
nario therein: with a longitudinal plasma density gradi- 
ent, the laser pulse can be forced to gradually increase its 
diameter, thus adiabatically turning from the pure bub- 
ble regime into the pure blow-out regime long before the 
laser is depleted. Fig. 1 shows a snapshot of a bubble 
around the transition from the bubble to the blow-out 
regime, where the electrons blown away from the laser 
and flowing around the bubble are so strongly deflected 
by the electric field of the captured bubble electons, that 
electrons are no more scattered into the bubble (which 
is not charge-neutralized yet). The number and abso- 
lute energy spread of the bubble electrons is thus frozen, 
but their energy is still increased due to the present bub- 
ble fields. The remaining energy of the laser allows to 
maintain the bubble structure for the remaining distance 
inside the capillary, where the laser is then guided along 
to overcome the Rayleigh limit. Only for this stage with 
increased laser beam size, the capillary-induced guiding is 
relevant. Note that this scheme is a two-stage approach, 
but within one and the same capillary, where the laser 
turns adiabatically from one stage into the other. It can 
be expected that the energy spread for 100 MeV is about 
1 % (as confirmed by measurements [2[), but 0.1% for 1 
GeV Normalized emittances are found both from 

simulations and experiments to range between 0.1 — 1 



mm-mrad. 



TABLE-TOP FREE-ELECTRON-LASERS 

An FEL requires an undulator which is an arrange- 
ment of magnets with an alternating transverse magnetic 
field. Electrons in an undulator are forced on a sinusoidal 
trajectory and can thus couple with a co-propagating 
radiation field. The induced energy modulation yields 
a current modulation from the dispersion of the undu- 
lator field. This modulation is called micro-bunching 
expressing the fact that the electrons are grouped into 
small bunches separated by a fixed distance. There- 
fore, electrons emit coherent radiation with a wavelength 
equal to the periodic length between the micro-bunches. 
In a Self- Amplification of Spontaneous Emission (SASE) 
FEL, there is no initial radiation field and the seed has 
to be built up by the spontaneous (incoherent) emission 

We present quantitative arguments, which are comple- 
mented by SASE FEL (GENESIS 1.3 0) simulations. 
The gain length, which is the e-folding length of the ex- 
ponential amplification of the radiation power, is 
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with undulator period X u and the basic scaling parame- 
ter p. This Pierce or FEL parameter describes the conver- 
sion efficiency from the electron beam power into the FEL 
radiation power and reads for the one-dimensional and 
ideal case (neglecting energy spread, emittance, diffrac- 
tion, and time-dependence) [3, [l6| 
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Here 7 = Eb eam /mc 2 is the electron beam energy, / 
the beam current, Ia = 17kA the Alfven-current, o~ x the 
beam size and A u = a u [Jo(() — Ji(()] (planar undula- 
tor), whereby a\ = K 2 /2, K the undulator parameter 
(K = 0.93 • X u [cm] ■ B [T] and B the magnetic field 
strength on the undulator axis), ( = a^/(2(l + a^)), and 
J are Bessel functions. In presence of energy spread, 
emittance, and diffraction, a correction factor A is intro- 
duced for the gain length [lj|: L gain = L gaintUeal (l + A). 
The saturation length L sat is the length along the undu- 
lator at which maximum micro-bunching is reached. The 
power of the FEL radiation at this point is the saturation 
power P sa t, while the shot noise power P n is the power 
of the spontaneous undulator radiation from which the 
SASE FEL starts. The coupling factor a = 1/9 describes 
the efficiency at which the noise power couples to the FEL 
gain process. The saturation length reads 
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where the saturation power scales as 
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The FEL- wavelength is, in case of a planar undulator, 
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Thus, for reaching a certain wavelength A, a shorter 
undulator period X u allows using less energetic electrons. 
In the following Tab. 1 and Fig. 2, we compare the 
FLASH VUV FEL (DESY) in the so-called femtosecond 
mode [13] with our corresponding table-top VUV FEL as 
well as a table-top XFEL operating with 1.2 GeV elec- 
trons, discussed further below. 



TABLE I: Parameters for the comparison between the DESY 
femtosecond-mode Flash-VW case and the corresponding 
table-top VUV FEL as well as a table-top X-ray FEL (rms 
values) . 



Parameter 


FLASH (fs) 


TT-VUV-FEL 


TT-XFEL 


current 


1.3 kA 


50 kA 


160 kA 


norm, emitt. 


6 mm- mrad 


1 mm- mrad 


1 mm- mrad 


beam size 


170 fim 


30 (an 


30 /im 


energy 


461.5 MeV 


150 MeV 


1.74 GeV 


energy spread 


0.04 % 


0.5 % 


0.1 % 


und. period 


27.3 mm 


5 mm 


5 mm 


wavelength 


30 nm 


32 nm 


0.25 nm 


Pierce par. 


0.002 


0.01 


0.0015 


sat. length 


19 m 


0.8 m 


5 m 


pulse length 


55 fs 


4 fs 


4 fs 


sat. power 


0.8 GW 


2.0 GW 


58 GW 



The importance of the ultra-high current in the table- 
top case is evident. The smaller undulator period allows 
a smaller beam energy, hence decreasing the gain and 
saturation lengths. The Pierce parameter gives the upper 
limit of the acceptable energy spread a^/j. Thus, for 
compensating the relatively large energy spread of laser- 
plasma accelerators and for maintaining a large output 
power, a beam current significantly above ~ 17 kA is 
mandatory for keeping p large and A small. 



Space charge effects 

Such ultra-high beam currents are subject to strong 
space-charge forces. After release into vacuum, the elec- 
tron bunch starts expanding, for which there are two 
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FIG. 2: Saturation lengths (Eq. [3} and (inset) degradation 
factor A as a function of electron beam current /: DESYs 
FLASH (dashed curve, triangle) and the table-top VUV sce- 
nario (solid line). The circles denote specific GENESIS runs. 



sources: (i) its own space-charge, driving a Coulomb- 
explosion, i.e. (transverse) space-charge expansion and 
(longitudinal) debunching, and (ii) its initial divergence, 
which drives a linear transverse expansion. Such ex- 
pansions transform potential energy into kinetic, hence 
changing the initial energy distribution and bunch form. 
First, we want to discuss an extreme case, that is, an 
upper limit of the charge that can be expected from a 
laser-plasma accelerator at lower energies. For this study 
we take a beam energy of 70 = 260, 1.25 nC charge, and 
a Gaussian bunch with sizes o~ x ,y,z — 1/im, correspond- 
ing to / = 150 kA, and an initial divergence of 6q = 1 
mrad. In the rest frame of the bunch, its length amounts 
to <t' z = er 2 • 70 = 260 /im. For such an aspect-ratio, the 
transverse Coulomb-explosion dominates over the longi- 
tudinal one (GPT [l8| simulations reveal that after a time 
period of 1 ps a' x — 73 fim, while a' z — 268 /im). Fig. 3 
shows the electron distribution in the bunch rest frame 
from the GPT simulations. 
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FIG. 3: (Color online) Projected spatial electron distribution 
within the bunch rest frame in the extreme case, at r = s 
(left) and r = 4 ps (right), from a GPT run. The color code 
indicates the electron energy (blue is initial kinetic energy, 
red is increased energy). The transverse expansion is clearly 
dominating. 



It can be seen that the transverse Coulomb explo- 
sion dominates the longitudinal one. We want to em- 
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phasize here that studies of space-charge effects of such 
ultra-dense relativistic bunches can be simulated most 
accurately only with codes (such as GPT and CSRtrack 
[ill) which utilize point-to-point interactions (a result 
also found in Ref. [20(). In contrast, other codes merely 
based upon Poisson-solvers cannot cope with large rela- 
tive particle motions within the bunch rest frame in case 
of large initial energy spreads and Coulomb explosion- 
induced motion. 

In comparison with the Coulomb-driven explosion, the 
linear divergence-driven transverse expansion is even 
larger: after a distance of 1 cm in the laboratory frame, 
the beam size is increased to 10 ^m, while transverse 
Coulomb-explosion yields 6 /jm only. Hence, even in the 
extreme case the transverse bunch expansion is domi- 
nated by the initial divergence and for lower charges this 
dominance increases. If 70 is the bunch energy just be- 
fore release into vacuum, 7' and j3' the energy and ve- 
locity of a test electron in the bunch rest frame, 7 and 
P « 1 its laboratory frame energy and velocity, then one 
can clearly distinguish between a transverse-dominated 
expansion and a longitudinal-dominated one: in the first 
case, where the electron moves in the rest frame with ve- 
locity ±/3 along transverse direction, 7 = 70 • 7', while 
in the latter case, where the electron moves with velocity 
±(3' in beam direction, 7 = 70 ■ 7' • (1 ± f3'). Hence, in 
a purely transverse expansion all electrons gain kinetic 
energy, while in a purely longitudinal expansion, some 
electrons lose energy. The inset of Fig. 4 shows the 
absolute kinetic electron energy distribution along the 
bunch at a distance of z = 3 cm after its release into 
vacuum (in case of 0.6 nC). This spectrum shows that 
the bunch under study undergoes mostly a transverse 
expansion. The maximum energy gain lies in the middle 
of the bunch, where the density is highest, i.e., where 
initially the highest potential energy was stored. That 
the space-charge driven transverse expansion is weaker 
than the divergence-driven one is also confirmed by the 
fact that (even in the extreme case) the divergence is 
increased only to 1.3 mrad. After passing a focusing sys- 
tem (a triplet of quadrupoles, positioned at z = 4 cm, 
having a length of about 10 cm), the Gaussian- like en- 
ergy distribution is transformed into a quasi-linear one, 
as depicted in Fig. 4. The reason is that in the result- 
ing quasi-collimated low-divergence beam faster electrons 
can catch up the slower ones. Note that after z = 3 cm 
(inset) the fastest electrons are symmetrically distributed 
around the middle of the bunch. These electrons start 
overtaking the slower ones, that is, they move forward 
within the bunch. If initially longitudinal debunching 
dominates, the fastest electrons are found in the head 
and the slowest ones in the tail of the bunch. 

A key parameter for SASE to occur is the energy 
spread 07/7, which must always be smaller than the 
Pierce parameter p. The relevant energy spread is taken 
over a bunch slice with the size of one cooperation 
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FIG. 4: Energy vs. position (in co-moving frame) along the 
bunch before entering the undulator (GPT simulation): the 
Gaussian-like distribution after a drift of z = 3 cm behind 
the gas-jet (inset), indicating a transverse-dominated bunch 
expansion, is stretched into a quasi- linear one behind the fo- 
cusing system. The vertical lines show the slice energy spread 
<t 7 and the horizontal lines mark the initial energy 70 = 260. 



length, i.e. the slippage length of the radiation along 
the bunch over one gain length. As seen in Fig. 4, 
a large fraction of the bunch fulfills 07/7 < p. Due 
to the slippage of the radiation relative to the bunch, 
the linear energy chirp must be compensated, because 
the energy factor 7 in Eq. |(SJ) increases along the 
bunch. However, in case of short-period undulators 
the gap is correspondingly small (typically one third 
of the undulator period) and such small gaps lead to 
strong wakefields also causing an energy variation along 
the bunch [24]]. In case of ultrashort bunches, the 
characteristic length of the resistive wake potential is 
longer than the bunch length and, hence, the energy 
change is found to be negative with a linear few-percent 
variation along the bunch. This in turn implies that 
the bunch as a whole is decelerated during the passage 
of the undulator, whereby the head loses less energy 
than the tail. It turns out that the linear energy chirp 
induced by space-charge corresponds well with the 
wakefield-induced slowing down of the entire electron 
bunch. Therefore, the radiation slipping in forward 
direction interacts with electrons of effectively constant 
energy, if by varying the gap the bunch energy loss is 
tuned with respect to the slippage. In other words, the 
space-charge effects and the wakefield effects cancel each 
other. We have found with simulations that for a first 
proof-of-principle of a table-top FEL the comparably 
low energy of 130 — 150 MeV is well suited to cope 
with the wakefields in the undulator. In this sense the 
table-top VUV case given in Tab. 1 can be regarded as 
an optimal test case. 

Coherent synchrotron radiation (CSR) [2l| as another 
source for increasing the (slice) energy spread within 
the undulator was found to have negligible impact. This 
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was shown with CSRtrack simulations which take into 
account that in our cases the bunches have a much larger 
transverse size (<J x ,y ~ 30 /nm) than length (a z ~ 1 /j,m). 



Table-top XFEL 

So far we have discussed a proof-of-principle scenario 
at relatively low electron energies, where space-charge 
effects play a dominant role leading to a linear energy 
chirp. The situation of the proposed table-top X-FEL 
(TT-XFEL) is different. For reaching a wavelength of 
A = 0.25 nm, an electron energy of f .74 GeV is needed 
in case of a period of X u = 5 mm. Space charge ef- 
fects are much weaker here. For this less demanding sit- 
uation we have confirmed with four different simulation 
codes (GPT, ASTRA fM], CSRtrack, HOMDYN fMj) 
that above f GeV, with the same parameters as in the 
extreme case given above, Coulomb-explosion leads to a 
projected energy chirp of below 0.3% and a bunch elon- 
gation with a factor below 1.1. Furthermore, CSR can be 
neglected in agreement with the ID theory. Experimen- 
tally the most demanding constraint is that the Pierce 
parameter is one order of magnitude smaller than in the 
test case at 130 MeV, hence the electron (slice) energy 
spread should be as small as 0.1%. This goal seems to be 
within reach considering the bubble-transformation sce- 
nario mentioned above. 

Without the effect of wakefields GENESIS simulations 
have shown that this TT-XFEL scenario with an undula- 
tor length of only 5 m yields 8T0 11 photons/bunch within 
~ 4 fs and 0.2% bandwidth, a divergence of 10 /xrad, and 
a beam size of 20 //m. However, the wakefields become 
the dominant degrading effect as the required undulator 
length is larger and the Pierce parameter reduced, hence 
also the tolerance with respect to variations in the elec- 
tron energy. But since there is no initial space-charge- 
induced energy chirp, one must find another method for 
compensating the wakeficld-induced energy variation. A 
suitable method for compensating the wakefields for the 
TT-XFEL would be tapering, i.e., varying the undula- 
tor period along the undulator. Due to the fact that the 
undulator parameter K is smaller than unity, tapering 
via K, i.e., by gap variation, could only be used as fine- 
tuning. Depending on the specific material properties of 
the undulator surface, our first wake calculations show 
that the bunch center loses over the entire 5 m undulator 
length about 10% of its initial energy. Due to the linear 
wake field variation along the bunch only a fraction of 
it will undergo SASE. A detailed study will be further 
investigated in a future paper. 



Experimental Status 

In order to demonstrate practical feasibility, we have 
built and tested a miniature focusing system and un- 
dulator consisting of permanent magnets. The focusing 
triplet consists of mini-quadrupoles with an aperture of 
just 5 mm, hence allowing for measured gradients of 530 
T/m. Their focusing strength was measured with a 600 
MeV electron beam (at Mainz Microtron facility MAMI, 
Mainz, Germany). A first test hybrid undulator with a 
period of only 5 mm and a peak magnetic field strength 
of ~ 1 T has been built and produced with a 855 MeV 
beam (also at MAMI) an undulator radiation spectrum 
as expected. These results will be published elsewhere. 



CONCLUSION 

We have shown by means of analytical estimates and 
SASE FEL simulations that laser-plasma accelerator- 
based FELs can be operated with only meter-scale un- 
dulators. The key parameter is the ultra-high electron 
peak current, which significantly reduces the gain length 
and increases the tolerance with respect to energy spread 
and emittance. The latter would also allow to increase 
the TT-XFEL photon energy into the medically relevant 
range of 20—50 keV, because the limiting quantum fluctu- 
ations in the spontaneous undulator radiation scale with 
7 4 [2^| and the required electron energy of the TT-XFEL 
is comparatively small. 
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